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Abstract 

We use analytic techniques and the dynamical mean field method to 
study the crossover from fermi liquid to polaron behavior in models of 
electrons interacting with dispersionless classical phonons. 
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I. INTRODUCTION 



In this paper we present the results of a study of the crossover from fermi hquid to 
polaron behavior in several related models of electrons interacting with dispersionless 
classical phonons. We use analytic techniques valid in weak and strong coupling limits, 
and we use the "dynamical mean field" method^ to obtain results at intermediate 
couplings. In a companion paper we present a detailed study of a more comphcated 
model believed to be relevant to the "colossal magnetoresistance" manganites. 

The electron-phonon problem has been extensively studied. The "polaron prob- 
lem" of a single electron coupled to a deformable medium has been understood in 
detail!. The problem of a degenerate fermi gas coupled to phonons has been solved 
perturbatively in a weak coupling limiti. There has also been fairly extensive work 
aimed at going beyond the original perturbative solutionis, but this work has been 
aimed mostly at understanding superconducting transition temperatures and charge 
density wave instabilities. The self-trapping or polaron physics has received less at- 
tention. 

There are several experimental motivations for our work. One is the "colossal 
magnetoresistance" materials Rei^x^xMnO^^, where Re is a rare earth such as lan- 
thanum and A is a divalent metal ion such as Ca or Sr. For 0.2 < x < 0.5 these are 
metals at low temperatures and "insulators" (in the sense that the resistivity is high, 
and rises as T is lowered) at high temperatures. The insulating behavior has been ar- 
gued to be due to a self trapping of carriers in local lattice distortionsEHHI. The physics 
of this material is complicated by the existence of an additional "double-exchange" 
carrier spin interaction, and will be discussed in a companion paper. La2-xSrxNiOi 
is another compound in which carrier localization by lattice distortions has been 
discussedS. Previous work on La2-xSrrcNi04 has focussed on density wave order at 
particular, commensurate, x but the effects discussed here are likely to be relevant 
also. We also mention the many semiconducting compounds which display polaron 
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effects when very lightly dopedll. At higher dopings, fermion degeneracy effects will 
become important, but the electron-phonon interaction will remain strong. 

We study models of electrons interacting with phonons. The electrons may or may 
not have spin and orbital degeneracy and are coupled to one or two independent local 
oscillators. The motivation for studying these different but closely related models is 
as follows. The interplay of polaron and fermi liquid physics is controlled by energy 
and entropy. The competing energies are the electron delocalization or kinetic energy 
Ekin and the lattice energy Eiatt gained by localizing an electron. The entropy of a 
localized electron depends upon the number of spin and orbital states available for it 
to localize into, while the entropy of the phonons depends on the phase space of the 
oscillators. Because the different models have different entropies, we expect them to 
behave differently at nonzero temperatures. 

Our study has several limitations. First, the numerical results are obtained via 
the dynamical mean field approximation. This is a controlled approximation which 
however becomes exact only in a limit in which the spatial dimension d oo0. 
Studies of other models have indicated that for = 3 the results are qualitatively 
correct and indeed in good numerical agreement with those obtained by other methods 
and by experiment00. The quantitative accuracy for quasi 2d materials such as 
La2-xSrxNi04^ is not clear. Further, the dynamical mean field approximation is 
essentially local. The only density wave instability which may be studied is the 
commensurate (vr, tt, tt, . . .) density wave, and the reliability of the approximation 
even for this case has not been established. We have not considered density wave 
solutions at all in this work. 

We have also assumed classical phonons. This approximation is known in the 
single-electron case not to affect the physics significantlyil (although it does lead to a 
few easily diagnosed low temperature pathologies), and it dramatically simplifies the 
computations. However, it does mean that we are unable to study superconductivity. 

A third limitation is that we study models in which the only interaction is 



the electron-phonon interaction. In particular, we omit the on-site "Hubbard- 
U" Coulomb interaction, which is undoubtedly important in at least some oxide 
materialsEi'0. While the approximations we have made limit the direct applicability 
of some of our results to experimental data on these materials, many of our results are 
experimentally relevant and we believe it is useful to have a detailed account of the 
behavior of a well-defined model in the literature. At various places in the text we will 
qualitatively discuss the effects of non-classical phonons and of the on-site Coulomb 
interactions, but a quantitative treatment will be left to a future publication. 

The remainder of this paper is organized as follows. Section II defines the models, 
the parameters, and the approximations. Section III discusses the numerical meth- 
ods used. Section IV presents the results of a detailed solution of the simplest model, 
a non-degenerate band of spinless electrons interacting with a single oscillator. We 
give phase diagrams, electron spectral functions, frequency, and temperature depen- 
dent conductivities and phonon probability distributions. Section V discusses the 
changes that occur if electron spin and orbital degeneracy, and multiple oscillators 
are included. Section VI is a conclusion. Appendix A clarifies the relation between 
the results presented here and those of the conventional Migdal-Eliashberg treatment 
of the electron-phonon problem. The other Appendices present details of various 
calculations. 



II. MODELS AND FORMALISM 
A. Hamiltonian 

In this section we define the models we study and present the approximations we 
use. We consider electrons coupled to phonons; the Hamiltonian, H, may be written 
as the sum of an electronic part, H^u a phonon part, Hph, and an interaction part 

Hei-ph as 
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H — Hel + Hph + Hel-ph- (1) 

For Hei we take electrons moving on a lattice. An electron on a given site i may be 
in one of Uorutai orbital states (we will consider Uorbitai — 1 and Uorutai — 2) and one 
of rispin spin states (we will consider Uspin = 1 and Ugpin = 2). The operator creating 
an electron of spin a in orbital state a on site i is and 

Hel^ - ^ ttjdtaadjbfi " X] dj^diaa- (2) 

ij iaa 

Here we have introduced the hopping matrix element and the chemical potential 

We model the phonons as localized classical oscillators. We write the displacement 
of the oscillator from some reference position as r. We consider two cases: "one- 
dimensional phonons", where r is a scalar quantity with — oo < r < oo, and "two 
dimensional phonons" where f is a two-component vector which we parametrize by a 
magnitude, r, and an angle, 0, via r = r{cos4>, sincf)). We assume that in the absence 
of the electron-phonon interaction the equilibrium position of the oscillator is r = 0, 
and that there is an elastic restoring force with force constant k. We measure r in 
units of lattice constant so r is dimensionless and k has the dimension of energy. We 
write: 

Hph = \kT''. (3) 

The form of the electron-phonon coupling depends on the physics. One possible 
couphng is of the phonon to the charge of the electron, i.e. 

Ht-X = 9Y.n{dtaadiaa - u). (4) 

iaa 

Here we have introduced the mean density n — -^{^iaadtaadiaa) and have defined 
r = to be the equilibrium phonon state for a uniform distribution of electrons. The 
electron-phonon coupling is g. If both the electrons and the phonons are degenerate, 
then one may have a coupling between the orbital state of the electron and the 
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direction of the local phonon. Such coupling may be expected to occur in materials 
such as Lai_xAxMn03 in which the Jahn- Teller effect is import antSEllll. In such 
materials there are in the absence of lattice distortions two degenerate electronic states 
on each site, transforming as a two dimensional representation of the cubic group. 
The degeneracy may be split by interaction with an appropriate phonon mode. The 
interaction expressing this physics is prescribed by group theory to be 

Hf_^, = gj:dl^^''d,,a-n. (5) 

iaba 

Here r = {tz,Tx) is a vector of Pauli matrices acting in the orbital subspace and 
r, = rcosd), = rsind). 



B. Parameters and Limits 

One important energy scale is the electron banding or kinetic energy. This depends 
upon the band filling n and the details of the hopping matrix elements t^j', but it is 
roughly of order t, where t is a typical value of tf-. 

Another important energy is that obtained by localizing an electron on a site. By 
minimizing Eqs. § or ^ and ^ one finds this energy to be of order g"^ /k. The ratio of 
these scales is the important dimensionless parameter A = g'^/kt; A < 1 corresponds 
to weak coupling and A > 1 to strong coupling. This A is shown in Appendix A to 
be related to the usual coupling constant of the Migdal-Eliashberg theory by a factor 

of ^spin^orb/ ^* 

A final important scale is the phonon frequency ujiy. In the weak electro n-phonon 
coupling limit it is known that for T > cu^i/S, phonons behave classically while for 
lower T quantum effects are important!. In the strong coupling limit, polarons are 
formed and large lattice distortions occur. From studies of a single electron coupled 
to a deformable mediumi one finds that quantum effects are important only at a 
very low scale uj£,/n, where n ^ 1 is roughly the number of quanta involved in the 
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formation of the polarona. Throughout this paper we assume that T is sufficiently 
high that quantum effects in the phonons may be neglected. 

C. Formalism 

Physical quantities may be calculated from the partition function, Z, which may 
be written as a functional integral over the electron and phonon fields. Because we 
have assumed classical phonons, we may solve the electron problem for given values 
of the phonon coordinate and then average over the phonon coordinates. Because 
the phonon is classical, we may rescale the phonon coordinate to r = \Jk/t r. The 
partition function becomes 



T-^phonon 



(6) 



2T^ 

Here T is the temperature, Vphonon is the phonon measure appropriate to the choice 
of phonon and G^ff is the electron Green function appropriate to Eq. |I] with fixed 
static Tj. Note also that r has been rescaled in such a way that in Hei-ph the coupling 
constant g = 

Starting from Eq. ^ one may make perturbative diagrammatic calculations using 
standard methods. In the weak coupling limit the starting point is rj = (corre- 
sponding to uncoupled electrons and phonons) and it may be seen from the first term 



in the exponent of Eq. y that fiuctuations in r are small, of order jT/t. This is the 



classical analogue of the Migdal's expansion parameter sjuio/t. In the strong cou- 
pling limit the starting point has some 7^ but fiuctuations are still of order \/T, 
so that at low T an analytic treatment is possible, if the ground-state configuration 
of the Ti can be determined. The connection to the conventional treatment of the 
electron-phonon problem is discussed in more detail in Appendix A. 

At intermediate coupling we have not succeeded in constructing the ground state 
about which to expand. We have therefore resorted to the "dynamical mean field" 
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approximationp in which it is assumed that the electron self energy is local (i.e. mo- 
mentum independent). In the electron-phonon problem in three spatial dimensions 
the self energy is known to have only a weak momentum dependence in the various 
solvable hmitJS, so the dynamical mean field approximation, which formally be- 
comes exact in a limit in which spatial dimensionality d oo, seems likely to be 
reliable. 

The electron Green function G in general depends upon orbital indices (a,b), spin 
indices {a, (3), momentum k and Matsubara frequency iun- The dynamical mean field 
Ansatz is 



Here e^^ is the Fourier transform of t^j*. All of the quantities in Eq. |^ are tensors 
in the direct product of spin and orbital space, and quantities without Roman (a, b) 
or Greek {a, (3) indices are proportional to the unit matrix in orbital and spin space 
respectively. 

Throughout this paper we shall assume that there is no long range order in ei- 
ther orbital space or spin space, so we may take ~ 6abSai3- Generalizing the 
formalism to include orbital ordering would be straightforward and quite interesting. 
The generalization to ferromagnetic spin ordering is given in a companion paper0. 
Because S is taken to be /c-independent all interaction effects are derivable from the 
local (/c- integrated) Green function, Gloc- Because we have assumed that there is 
no long-range order in orbital space, G^oc ~ ^ah- We may therefore write 



The k integral may be simplified by introducing at each fc-point the ah space rotation 
Rk which diagonahzes G""^, writing G°'^{k) = RkG DiagRk^ , and exploiting the cyclic 
invariance of the trace. The k integrals of the two components of Goiag must be 
identical, so we may finally write 




(7) 
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GLOC{l^n) = / -. ^^7^ — ^ (9) 



with V{€k) the density of states. 

Because all interactions are local, Gloc niay be derived from a partition function 
identical in form to Eq. ^ but with spatial index i suppressed and the quantity 
replaced by 

G MF ~ d^a^aa + Hf.l-ph (10) 

at fixed phonon coordinate. The free energy corresponding to a given value of phonon 
coordinate may be found as usual from Gmf, and from this one may construct the 
phonon probability distribution 



P(r) = — exp 

Zloc 



^r^ + Y.Tr In Gljp 



where the local-approximation partition function Z^oc is defined by 

VphononPir) (12) 



^LOC = 

Here a is a frequency-dependent effective field (analogous to the magnetization in the 
usual Weiss mean field method) determined by the condition 

d In Zloc ..o^ 
GLOcyi^n) = -^—r- — ^- (13) 

This condition may be implemented by observing that 

T,{iujn) = a{iixJn) - Gloc- (14) 

The precise form of the resulting equation depends on the density of states. We shall 
consider two forms, the Lorentzian 



T-^ Lorentziani^^k) ~n \ TTT (1^) 

+ 



and the semicircular. 



Vsemi{ek) = ^^t^ -ell {27^ e). (16) 

The Lorentzian density of states has the advantage that the self consistency equation 
may be analytically solved to yield 

aLORii^n) =iuJn + fJ' + itsgnUn- (17) 

It has however two unphysical features: 'Diori^k) is unbounded and 
/^o^ dekekT>Lor{^k) = C)0 so the kinetic energy is infinite for finite fi. 

For the semicircular density of states the integral may also be performed ana- 
lytically, leading to 

asemiiii^n) = iuJn + - G LOcii^n) ■ (18) 

This equation must be solved numerically, except in simple limits. 

The semicircular density of states corresponds to a Bethe lattice and the dynamical 



mean field theory together with the self-consistency condition Eq. |18| can be derived 
in this case by a more physically transparent argument. We take for concreteness 
the case of orbitally degenerate electrons with a Jahn- Teller coupling to orbitally 
degenerate classical phonons. On a Bethe lattice the local Green function Gu on a 
site i with phonon amplitude is: 

Gt^\^^n,n) = iu^ + ^i + gT-r^^-Y.npf^{iu:n,r,)t% (19) 

jcd 

One may diagonalize the Jahn- Teller coupling by a local rotation Ri in orbital 
space, with Ri = cos(f)i/2 + sin(j)i/2 r^, and one may absorb Ri into the hopping 
matrix tij = RitijRj which thus acquires a dependence on the relative orientation of 
the phonon distortions. The mean field theory for the disordered phonon state follows 
(in the limit of large coordination number) from the assumption that the 0j's and rj 
are statistically independent and the sum on j can be replaced by an average over 
0's and r's, the latter with a non-trivial distribution P(r) defined by Eq. 0. The 
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hopping self-energy then reduces to S'^'^Gloc with Gloc = 1/2 Tr < Ga and the 
self-consistent equation for Gloc is obtained by averaging Ga given by Eq. |19| over 
Tj. Eq. |1^ may be helpful as a starting point for mean field theories with non-trivial 
spatial correlations of the phonons. 

We will be interested in three physical quantities: the momentum-integrated elec- 
tron spectral function A{uj), the phonon probability distribution P{r), and the con- 
ductivity a. The spectral function is 

TT 

In the dynamical mean field method the conductivity (Tab = o'MF^ab with 

aMF{i^) = J depV{ep)Tj2G{p,iuj)G{p,iuj + in). (21) 



In particular, the dc conductivity is 

adc{T) = eH^Tr J ^P(ep)[/m G{ep,uj)]y4Tcosh{uj/2Ty (22) 

with G(ep, u)-^ = u} + ^- (ef ) - E(cj). 

As T — > 0, two possibilities arise: either E"(u; = 0) or it tends to a finite 
non-zero value. In the former case, 

a,.^et n^^^n,^bV{e ) j ^J^^^^^Ji^^ (23) 
with e* satisfying e* = /i — S'(0). In the latter case 

S"(0) 



TT 



(ep + /i-S'(0))2 + S"(0)^ 



(24) 



III. NUMERICAL METHODS 

In this section we outline the numerical methods used to solve Eq. and compute 
physical quantities. To aid in the discussion we write this equation explicitly for the 
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simplest case of nondegenerate electrons coupled to nondegenerate phonons. We have, 
setting t = 1, 

nOO P{r) 

aiuj) = 100 + 11- / dr , , — (25) 
J-oo a{uj) +gr ^ ^ 

with 



P(r) = — exp 

^LOC 



„2 



2T ^ 



a„ + gr 



(26) 



and Zi^oc given by Eq. |T2|. Here = Q.{%u)n) with a;^ a fermion Matsubara frequency 
LOn = (2?7, + l)7rT. 

For most physical quantities we require a{uj) for real frequency. It is possible to 
express P(r) in terms of a{uj), but this involves performing numerically an integral 
over a continuous frequency. We have found it more convenient to find P(r) by solving 
Eq. ^ on the Matsubara points and then to use this P{r) to solve for a{uj) on the 
real axis. In what follows we first discuss the solution on the Matsubara axis and 
then mention the additional issues arising for the real axis case. 

We solve Eq. ^ on the Matsubara axis by direct iteration, so that a„ at step 
+ 1 is determined by evaluating the right-hand side of Eq. ^ using a„ from step 
N. If a solution at the same coupling and a nearby temperature is available we use 
this as the starting point; if not we use the g = solution. Convergence is usually 
rapid; typically 10 iterations are required for convergence from the g = solution to 
a solution with RMS error of 10~^ (averaged over all retained Matsubara frequencies) 
and using the solution from a nearby T saves about 2 iterations. However away from 
half filling, in the strong coupling limit, the straightforward iteration procedure does 
not converge, but instead goes to a two-step limit cycle in which evaluating the right 
hand side of Eq. |2^ using the values a^^ yields a different set of values a^^^ which, 
when put back in Eq. generate again This limit cycle may be avoided by 
making the iteration proceed in smaller steps by defining a„ at step + 1 to be an 



appropriate linear combination of the right-hand side of Eq. ^ and the solution at 
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step N; however as T gets lower or coupling gets stronger the required admixture of 
the newly computed a gets smaller, and at some point the computations become too 
time consuming. 

Away from half filling it is necessary at each temperature and coupling to find 
the chemical potential fi{T,g) which gives the desired n. We do this by computing 
at several values of fi and interpolating. At strong coupling and low T, n is a very 
sensitive function of /i, and it is necessary to choose at least three /i- values all of 
which yield n's which are within 10 percent of the desired value and to use a cubic 
interpolation. In the calculations presented n varies by less than 6 x 10^^ over the 
range of T and g considered. We found it convenient to evaluate n via 

n = hm^^oT J2 Tr GLod^^n) e'^^-^ (27) 

n 

Gloc may be expressed in terms of a via the mean field equations. 

The next issue is the number of Matsubara points needed to compute P(r) accu- 
rately. We note from Eq. |2^ that at large u;„ 

an = iuJn + fi-- \ fdr P{r){fi + gr) + (^) (28) 

Analogous but slightly different formulas apply in the case of Jahn- Teller coupling. 
We determine the maximum retained Matsubara frequency rimax by requiring that 
the 0{l/ujfi) terms are less than 10~^ and in the expression for P(r) we evaluate the 
terms with \n\ > rimax analytically, using the a„ given in Eq. |2^; the errors are of 
order n~^^. Equation ^ is also useful for improving the accuracy of the computation 
of the density n from Eq. ^ away from half-filling. To achieve the stated accuracy 
we found it necessary to choose Umax so that Un^ax = '^'^Tirimax + 1/2) is about five 
times the full width of the real-axis electron spectral function. The 1/T growth of 
^max as T — > limits our ability to calculate at very low temperatures. For example 
at ^ = 1.29 and T = 0.02, n^ax = 250. 

The remaining numerical issue for the Matsubara-axis solution is the integral over 
the phonon coordinate, which is done using the numerical quadrature routine "quad" 
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from the "Port" library. Gaussian quadrature routines are faster at weak coupling, 
but fail at strong coupling and low T where P(r) has a large narrow peak at an 
a-priori unknown location far from the origin. 

We now turn to the real axis calculations, which we also do by direct iteration of 
Eq. ^ with analytically continued frequency uj + ie, using P(r) which has already 
been computed. The limit cycles which plagued the Matsubara calculations do not 
occur, but other numerical problems arise. Because P(r) is held fixed, different 
frequencies decouple and the mean field equation does not come from minimizing a 
free energy. Nearby frequencies are also sometimes found to converge at different 
rates. In parameter regimes in which the T = spectral function has a gap, the 
imaginary a' {uj) part of a{Lo) is small and very temperature dependent at low T for 
frequencies in the gap. Our procedure fails to find a non-zero value of a" if this 
quantity becomes less than 10~^, and we are not sure of the accuracy for a" < 10~^. 
Physical properties in such regimes depend crucially on the tailing of the spectral 
function into the gap, so our ability to calculate the low T properties of regimes with 
gaps is limited. 

To calculate conductivities we solve the real axis mean field equations for a discrete 
set of frequencies spaced by 0.01. Then at these frequencies we evaluate 

e{uj) = J depV{ep)[ImG{ep,uj]'^ (29) 

using the Port routine "quad" and then obtain 

/UJ 
du 6{uj) cosh""^— (30) 

using the quad routine and evaluating 9{uj) by spline interpolation from the known 
values. In the computation we set the factors e^/vr = 1. 

We note finally that especially in the crossover region where there is a frozen-in 
lattice distortion which is not large enough to open a gap in the spectral function, 
physical properties are quite sensitive to the precise value of the coupling constant. 
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A change in g by 2% can lead to more than factor of three changes in e.g. the 
low T conductivity; also small numerical errors in the computation of the fermion 
contribution to the action can mimic a small change in g. An example can be found 
in the curve labelled g = 1.60 in Fig. 2. The analytically calculated T = resistivity 
for this g is 3.75. The actual numerically calculated low T spectral function is that 
expected for g = 1.57; the corresponding analytical T = resistivity is found to be 
1.004, in good agreement with the numerical calculation. 

IV. NON-DEGENERATE SPINLESS ELECTRONS 
A. Introduction 

In this section we present results of a detailed study of a model of spinless, orbitally 
non-degenerate fermions coupled to classical non-degenerate oscillators. This section 
is organized as follows. In this introduction, we outline the qualitative physics. In 
part 4-B we solve the original model, Eqs. [l|, 0, |^, ^ in simple limits, obtaining results 
which clarify the interpretation of the d = oo results. In part 4-C we solve the d = oo 
model with the Lorentzian density of states, (Eqs. |T^, |T7p, and in part 4-D the d = oo 
model with semicircular density of states (Eqs. [T2|, |TB|). 

In the model defined by Eqs. 1-4 we may distinguish three phases according to 
the T — >■ limits of the oscillator coordinate r and the electron spectral function. 
In weak coupling, limT-^o{r) = and the spectral function takes the non- interacting 
form. The low T resistivity is then linear in T, and extrapolates to as T 0, 
being due to the thermal fluctuations of r about r = 0. In intermediate coupling, 
limT^Q{r) 7^ 0, implying that some frozen-in lattice distortions exist. However, the 
amplitude of these distortions is too small to localize states near the fermi level. The 
low-T resistivity is still linear in T, but extrapolates to a non-zero value at T = 0. 
The spectral function is perturbed from its non-interacting form. In strong coupling 
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the amplitude of the frozen-in lattice distortions is large enough to localize electrons 
at the Fermi level, leading to a gap in the T = spectral function and an activated 
resistivity. 

All three regimes are found in the calculations discussed below. None of the 
regimes involves long rang order- the localization should be thought of as being due 
to polaron formation. At commensurate filling or in low dimensions, effects of long- 
range order (i.e. charge density wave formation) will be important but at general 
incommensurate fillings these effects may be neglected. 

Quantum effects will also be important: at sufficiently low T and in the absence 
of long range order, quantum tunnelling will restore translational invariance. In 
the strong coupling limit the quantum effects may be thought of as leading to the 
formation of a polaron band, with a band width exponentially small in A; at scales 
greater than the polaron band width the results presented here will apply. The effect 
of quantum fluctuations on the intermediate couphng regime is less clear, but must 
be left for future work. 



B. Full Model, Solvable Limits 



The crucial issue is performing the integral over the phonon coordinates in Eq. ^ 
At low T the integral is sharply peaked about particular values {r*} and a steepest 
descents approximation is possible. Two solvable cases occur: 

1. Weak coupling. For A < Ac and T ^ (Ac will be estimated below) r* = and 
one may expand the logarithm in Eq. ^ obtaining 

(1 - XBik)) nr^k' 



Z oc exp 
Here Gq^ = iujn 



Tr In Go ^ 



Vr-k exp 
ek + and 

J (2ttY e„ - ep+k 



2T 



(31) 



(32) 
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It is evident from Eq. |3T| that the expansion in powers of r about rj = is controlled 
at low T if A < Ac with 



= maxk B{k). (33) 

In particular, as T — > at fixed A, r^^ ~ \/T is small, so the leading order 
perturbation expression for the fermion self energy, S, applieJi@. This expression is: 

E(p,iuj,T) = XT / - — ^ — -— . 34 

' ^ J (27r)3 1 - Afi(A;) icj - ep+fc ^ ' 

For small ujn, we may restrict the k integral to wavevectors such that p' = p + k is on 
the Fermi surface. The leading contribution to the u dependence of S is then 

S(p, iu, T) = inTX sign iun 

I ^V{e,,)/[l-XB{p-p')] (35) 

and if also A <^ Ac 

T,{p,iu!,T) = mTXV^ep) sgnuin (36) 

with f){ep) the angle averaged density of states at the fermi surface. 

Equation |36| is shown in Appendix A to be precisely the usual Migdal expression 
for the self energy. If (Ac — A)/Ac is not too small, the quantity AT / P(ep/)/[l — 
XB{p — p')] depends only weakly on p; if it is averaged over p or the Fermi surface it 
becomes the classical limit of the usual "a^F" function. Here T /t plays the role of the 
usual Migdal parameter un/Ep. As A — Ac, the phonon propagator becomes large 
and also acquires a singular momentum dependence because the phonon propagator 
diverges at a particular wavevector, so one may expect the Migdal approximation to 
fail. A discussion of corrections to the Migdal approximation is beyond the scope of 
this paper. Note also that A < Ac is the condition for the linear stability of the = 
state. A T = first order transition to a state of rj 7^ at a A* < Ac is possible and 
indeed will be demonstrated in the dynamical mean field approximation. 
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The resistivity may also be computed. If A ~ Ac, the momentum dependence of 
the phonons may be important and vertex corrections must be considered. If A ^ Ac 
the momentum dependence of the phonons, and therefore the vertex corrections, may 
be neglected and the conductivity on the Matsubara axis is 

■ G{p, iun + iQn)G{p, iun). (37) 

Substituting Eq. ^ into Eq. performing the integrals and analytically con- 
tinuing gives 

In a simple tight binding model, Vx{p) = 2tsinpx and evaluating Eq. |3^ we have 

TT AT 

2. Full Model, Strong Coupling. In the limit A ^ Ac an analytic solution is possible 
because we may treat the electron kinetic energy as a perturbation. If it is neglected 
then to solve Eq. 1-4 we note that if the mean density is n there are n sites with one 
electron, and 1 — n sites with no electron. On the n occupied sites there is a lattice 
distortion of amplitude 

r:cc = -i9/k)il-n) (40) 
and on the 1 — n unoccupied sites 

r*unocc = i9/k)n. (41) 
The energy gained is found from Eq |l] to be 

Egained = -\t\{{l - uf + n'). (42) 

Perturbative corrections (in powers of the kinetic energy) may be calculated. They 
involve processes in which electrons make virtual hops to unoccupied neighboring 



sites, and lead to a repulsive short ranged interaction of order t/\ between polarons. 
This interaction will lead to polaron ordering, at a scale generically much lower than 
t/X. 

In this limit transport is thermally activated; the prefactor is t and the gap is 

■^gained- 

C. Dynamical Mean Field Method: Lorentzian Density of States 

In this subsection we discuss the solution of the dynamical mean field equations for 
spinless fermions with Lorentzian density of states coupled to a single scalar oscillator. 
We begin with the T ^ limit, where one may use a steepest descents approximation. 
The details are given in Appendix B. We find that for A < Ac(n), r* = lirriT-^oi^) — 0> 
while for A > Ac there are two extrema, at r* < and > 0. For n = 1/2, rjj' = — 
the transition at Ac is second order; and r* grows smoothly with A. For n ^ 1/2, 
rl 7^ —r^ and the transition is first order. 

The physical content of the r* 7^ states is most clearly seen from the electron 
Green function which at T = is 

GLOc{^u^) = . ^ ^.f — ^ + . ^ ^.f — ^ (43) 

too + jj, + ttsgnuj + grl too + jj, + itsgnuj + gr^ 

with Zi = Ji^l «| and Z2 = , J,'^! , Thus if r* 7^ the spectral function consists of 

''ll + l''2l F1I + F2 

two resonances. In the large A limit, these are well separated. One, of weight n, is at 
the negative energy fi — g{l — n) and corresponds to the sites on which the electron- 
phonon interaction has localized an electron. The other, of weight 1 — ra, is at the 
positive energy ^ + gn and corresponds to the empty sites. The splitting, AE oc A as 
expected, while the energy gained Egamed (relative to the undistorted state) is 

Egained = -^At[(l - flf + U^] (44) 

consistent with Eq. The Lorentzian density of states is pathological in that, as 
may be seen from Eq. at any value of g the T = spectral function is non-zero at 
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uo = jjL] thus the ground state is always metalhc. This pathology may be traced back 
to the non-integrability of the first moment of density of states, which implies that 
some carriers have arbitrarily large velocities and cannot be localized by any lattice 
distortion. 



D. Dynamical Mean Field Method: Semicircular Density of States 

In this subsection we present results obtained using the dynamical mean field 



method with semicircular density of states. Our analysis is based on Eqs. |12|,|T8|. The 
technical differences from the previously discussed Lorentzian case are that the mean 
field function a{uj) satisfies a non-trivial self-consistency equation, and the moments 
of the density of states are finite, so the spectral function may develop a gap. 

As in the previous subsection, one may use steepest descents arguments for T 
near zero. The details of the analysis are given in Appendix C. For g < g* the energy 
minimum is at r* = 0; for g > g* there are two minima, at r = ri > and r = r2 < 0. 
For n = 1/2 the transition is second order and occurs at g* = g^ = v^37r/2 = 1.535; 
for n 7^ 1/2 the transition is first order and occurs at a (7* > gc- 

As in the previous section, the physical content of the r* 7^ solutions is the 
existence of frozen-in lattice distortions. The effect of these may be seen from the 
local Green function which, at T = is 

Gloc = { , + ) /(^i + ^2) (45) 

\a(u;) ^ grx a{uj) + gr2 ) 

Here Zi > and Z2 > are the relative weights of the free energy minima corre- 
sponding to r = ri and r = r2 respectively. 

Figure 1 shows A{uj) = —IitlGloc for n = 1/2 and different values of ri = r2 = r. 
Our interpretation of these results is that for < r < 1/2 the potential fluctuations 
due to the frozen-in phonons are very weak, and merely localize a few states at the 
band tails. This accounts for the slight broadening of A and the slight decrease near 
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Lo — —/I. For 1 > r > 1/2 the potential is strong enough to create a minimum in A 
at a; = — for r > 1 it is strong enough to locahze all of the electrons and create 
a gap in the spectral function. Note also that as shown in the inset to Fig. 1, r 
is rapid function of A in the region r < 1. This rapid growth of r means that the 
transition is very nearly first order even for n = 1/2. Figure 2 shows the temperature 
dependence of the resistivity for the same coupling constants used to construct fig. 1. 
The two lowest resistivities are calculated for weak couplings g < Qc and g — Qc, for 
which limT^or*{T) — (so the T — spectral functions are identical). For g < gc, 
p{T) at low T is linear in T with g(-dependent prefactor and limT^opiT) = 0. For 
g — gc the quadratic term in the phonon energy vanishes and p ~ T^/^ at low T. 
For g slightly greater than g^ lirriT^o p{T) ^ but is finite. For g large enough 
that the spectral function has a minimum at a; = 0, the resistivity initially drops as 
T is raised. For g large enough that the T — spectral function has a gap, p rises 
rapidly, and ultimately exponentially as T — > 0. The temperature at which p begins 
to rise rapidly is somewhat less than the T = gap. Note that for \/37r/2 < g < 1.63 
limT-*Qp{T) — p(0) is neither zero nor infinite. 

Some insight into the process of gap formation can be gained from Figs. 3 and 4, 
which show the temperature evolution of the spectral function and phonon probabil- 
ity distribution for coupling g — 1.69. One sees from fig. 3 that the rise in resistivity 
as T is decreased begins at the T — 0.07 at which the spectral function first develops 
a minimum. A weak maximum in P{x) is visible in fig. 4 even at higher temper- 
atures, but the maximum becomes an obvious feature at the temperature at which 
the resistivity begins to turn up. The initial rapid rise in the resistivity is associated 
with the development of the minimum in the spectral function; the asymptotic low T 
behavior p ~ e^^'^ occurs at or below the lowest temperatures available numerically. 
The numerical calculations in this section have all been performed at n = 1/2, where 
particle- hole symmetry implies p = 0, simplifying the computation. Results for the 
particle-hole asymmetric case in other models will be presented below. In our com- 
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putations to date we have found that n ^ 1/2 has quahtatively similar behavior to 
n = 1/2 with metaUic and insulating regimes according to whether or not the T = 
spectral function has a gap. As shown by the explicit calculations of Appendix C, 
somewhat stronger couplings are required to obtain insulating behavior at 7^ 1 than 
at n = 1 In the numerical calculations at T > 0, we have seen no evidence of the first 
order transitions found at T = 0. 

E. Effect of Electron and Phonon Degeneracy 

N on- degenerate phonons. The effect of electron spin and orbital degeneracy on models 
with scalar phonons is straightforward, if there is no long range order. One multiplies 
the quantity Tr ln[G~^j:] in Eq. ^jby the factor Uspinnorb- This factor may be absorbed 
into T and g by defining T = Tuspmnorbi and g = g / y/n^^^i^jT^ and the model reduces 
to the previous case, apart from small differences coming from the T-dependence of 
the fermion action. This is shown in Fig. 6, which depicts p(T) for a model with a 
two- fold electron degeneracy {rispinnorb = 2). The couplings are chosen to differ by a 
factor of \/2 from those used to construct Fig. 3, so as to produce the same T = 
behavior. The resistivity is smaller by a factor of two, because it is proportional to 
1 / Uorhnspin- If the temperature axis is rescaled by a factor of two, the various curves lie 
almost on top of each other, with agreement being best for low T and weak coupling. 
For example, the temperatures at which the resistivity turns up for g = 1.83 and 1.69 
in the non degenerate case are almost exactly one half of the temperatures at which 
p turns up for g = 1.30 and 1.20. Some deviations from this scaling become apparent 
at higher T; these are due to the T-dependence of the fermion action (i.e. to the 
electron entropy). For example, the scaled weak coupling resistivities differ by an 
amount proportional to T^. Similarly, the upturn for the strongest coupling occurs 
at a slightly lower rescaled temperature in the degenerate case. Again for g ^ gc the 
model is insulating as T — 0. The approach to the T — limit for this model is 
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shown in the middle panel of Fig. 5. It is seen that in this model our calculations 
can access more of the low-T limit because of the rescaling of the temperature. 
Degenerate Phonons Jahn-Teller Coupled to Orhitally Degenerate Electrons. The 
theory is based on Eq. 5; significant differences from the previous cases occur because 
the electron spectral function may have a more complicated structure and the phonon 
measure suppresses the probability of small amplitude lattice distortions and increases 
the probability of large amplitude distortions. 

We begin with the effect of the phonon measure. For this discussion we restrict 
attention to n = 1, so the only difference in the model is the factor of rdr in the 
phonon measure. Figure 7 shows the resistivity of this model for n = 1 at the 
same coupling constants used to construct Fig. 6. One sees immediately that all 
resistivities are larger in the Jahn-Teller case; also for parameters such that there is 
a T = gap, up-turns occur at higher temperatures. The origin of these differences 
may be seen most easily in the phonon probability distribution shown in Fig. 

8 for both degenerate and non degenerate phonons using in both cases degenerate 
electrons. We have chosen g — — 1.085 (the critical value at which the T — 
lattice distortion vanishes) and T — .15. The results are representative of all g and 
T. One sees immediately that the mean square value of the lattice displacement is 
larger in the Jahn-Teller case than in the non-degenerate case and more importantly 
the small r fluctuations are suppressed. The larger mean square displacement means 
more scattering and hence more resistivity. The suppression of small r fluctuations 
means that P(r) is reasonably sharply peaked at a non-zero value, r = Tpeak- The 
temperature dependence of Vpeak in the degenerate-phonon case is determined by the 
couphng; for g < gc TpeakiT) ~ T^/^ for g = gc TpeakiT) ~ T^/-^ and for g > gc Vp^ak 
tends to a constant. 

For non-degenerate phonons and 7 = 0, P{r) has a maximum at r = if < g^c 
and maxima aX r ^ {) ii g ^ g^. Even for g < g^ P{r) aX T ^ {) may have a weak 
maximum at r 7^ 0, but this does not seem to have any consequences for physical 
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properties. In the degenerate-phonon case the peak structure of P(r) actually can 
lead to a minimum in the spectral function at zero frequency if grp^ak exceeds a critical 
value somewhat larger than the 0.5 value found in Appendix C. Of course if (7 < Qc, 
this minimum will vanish as T is decreased. This behavior is illustrated in Fig. 9, 
which shows the temperature evolution of the spectral functions for the two models. 
Only the degenerate phonon case has a minimum at u = 0; the minimum vanishes 
below T ~ 0.045, at which temperature grpeak ~ 0.5. This may be understood by 
reference to the T = limit, in which P{r) may be approximated by a delta function 
at a coupling dependent value tq. If gro > 0.5 the T = spectral function develops 
a minimum at = 0; if gr^ > 1 it develops a gap. Similarly, at T > the peak in 
the degenerate-phonon P(r) leads to a minimum in A{uj) if this peak occurs at an 
Tpeak such that grpeak > 0.5. This tendency to open or increase a gap in the electron 
spectral function as T is raised acts to increase the resistivity of the degenerate model 
above that of the non-degenerate one. 

Finally, the approach to the T ^ limit for the insulating regime is shown in 
the upper panel of Fig. 5. It is interesting that the result approaches the T 
limit more smoothly than in the two previous cases, and also that Tlnp is larger at 
T > than at T = 0, unlike the other two cases. We speculate that the origin of this 
difference is the larger rms value of r, which would imply an effectively larger gap. 

Away from n = 1 a more important difference in physics occurs, which is most 
easily seen for spinless electrons in the very strong coupling limit. Suppose n < 1 
and neglect hopping completely. Then by analogy with the steps leading to Eq. ^ 
one expects n sites occupied by a single electron and 1-n unoccupied. The spectral 
function associated to the singly occupied sites has two peaks, one below the chemical 
potential corresponding to the occupied orbital and one above corresponding to the 
unoccupied orbital. The spectral function associated with the unoccupied sites has 
only one peak, because in the absence of any electrons there is no lattice distortion and 
therefore no splitting of the orbitals. Consider now the first perturbative correction 
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due to the hopping. This will lead to a charge density 5n ~ n/A on previously unoc- 
cupied sites. If (as assumed in our application of the dynamical mean field method) 
there are no intersite correlations in orbital occupancy or phonon fluctuations, this 
extra charge will be randomly distributed over the two on-site orbitals. Thus in the 
presence of a small amplitude phonon distortion 5r the charge density will lead to an 
energy gain of at most 5n(5r)^, which will be too weak to compete with the phonon 
stiffness in the strong couphng hmit. The three peaked spectral function is therefore 
stable in the strong coupling limit. Appendices D and E derive these results from 
an asymptotic analysis of the dynamical mean field equations. It is clear from the 
above arguments that the structure of the mid-gap states may be affected by intersite 
correlations, which could lead to a charge fluctuation favoring one particular orbital, 
and therefore to an energy gain proportional to 6r rather than 6r'^. Incorporating 
intersite correlations is an important open problem. 

Of course, as A is decreased, the possible energy gain increases and at some point 
the mid-gap states split and the two pieces join the upper and lower bands. The 
resulting intermediate coupling regime corresponds to a uniform electron density and 
a uniform nonzero lattice distortion. The electron spectral function consists of two 
bands, with the fermi level in the lower one (for n < 1). The rigid band like state 
gains kinetic energy ~ {1 — n) relative to the three band state, because one band is 
not fully flUed, but is shown in Appendices D and E to be more costly in Jahn-Teller 
energy because the lattice distortion is not as large. 

As the coupling is increased at T = three phases are in principle possible: a 
weak coupling phase with a single-peaked spectral function, a "rigid band" phase 
with a two peaked speactral function and the fermi level in the lower band, and 
an insulating three peaked spectral function. Determining the couplings at which 
the system goes from one phase to another in general requires numerical solution of 
the coupled mean field equations, along with comparison of the energies of different 
locally stable solutions. However, in the limit {1 — n) the equations simplify and 
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a straightforward physical argument may be used. Consider the n = 1 problem in the 
insulating regime in which the spectral function has two peaks separated by a gap 
and calculate the electron removal spectrum allowing the on-site lattice distortion to 
adjust. (The u < part of the spectral function previously computed for this problem 
represents the lowest energy way to remove an electron from the system at fixed lattice 
distortion.) The energy may be computed from the quantity G1^{uj,r) defined in Eq. 



19| , with the G on the right hand side of the equation given by the n = 1 solution 
of the dynamical mean field equations. Performing this computation and optimizing 
over r shows that there are two extrema-one at r = and one at r = r* (with r* 
the optimum r found for n = 1). If g > 1.308... (i.e. r* > 1.527... ) then the energy 
of the r = extremum is lower, showing that the three peaked spectral function is 
favored, while for g < 1.308... the "rigid band" solution is favored. (Note that at 
g = 1.308 the two peaks of the rigid band spectral function are spearated by a gap). 
For this reason we believe that for (1 — ?^) ^ 1 there is as g is increased a second order 
transition to the rigid band two peaked spectral function. As g is increased further, 
the peaks separate and eventually a gap develops between them (although the ground 
state remains conducting). This is followed by a first order transition at a larger g to 
an insulating state with a three peaked spectral function. This sequence of spectral 
functions is shown in Fig. 10. On the other hand, for n ^ the results of Appendix 
E show that the weak coupling one=peak spectral function proceeds directly to the 
three peaked one via a first order transition. 

One important consequence of the central feature appearing in the three-peak 
regime is that it acts to fill in the gap created by the lattice distortion, so that one 
must go to stronger couplings and lower temperatures to see insulating behavior in 
the Jahn- Teller case than in the non-degenerate case. This may be seen directly from 
the strong coupling calculations of Appendices C and E. In the non-degenerate case 
the spectral function has two features, at ±(7^ and physical quantities are controlled 
by the energy gap between them, which is of about the same size. In the Jahn- Teller 
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case the highest and lowest of the three peaks are still separated by an energy of 2g^, 
but physical processes are controlled by the gap separating the lowest and middle 
peaks, which is much less. This may also be seen in numerical calculations of the 
low T spectral functions and phonon probability distributions. The upper panel of 
Figure 11 compares the spectral functions for the non-degenerate and Jahn-Teller 
cases at T = 0.03, g = 1.58 and n = 0.75. Although the true g > 1 strong coupling 
limit has not been reached, the two peak and three peaked structures expected for 
the non-degenerate and Jahn-Teller cases are clearly visible and the gap in the Jahn- 
Teller case is much less. Similarly, the lower panel compares the spectral functions 
for n = 1 and n — 0.75, for g — 1.58, T — 0.03 and Jahn-Teller couphng. Again 
one sees the effect of the mid-gap states. Figure 12 similarly compares the phonon 
probability distributions in the two cases. Figure 13 compares the resistivities of 
the two models for the weaker coupling g — 1.29; the non-degenerate phonon model 
displays insulating behavior below T ~ 0.13, while the Jahn-Teller model remains 
metallic. The resistivities cross at higher T because of the larger fluctuations of (r^) 
in the Jahn-Teller case. 



V. CONCLUSIONS 

We have used the dynamical mean field method to study the crossover from 
Fermi hquid to polaron behavior in models of electrons coupled to locahzed classical 
phonons. The models we studied involved electrons with and without spin and orbital 
degeneracy coupled to degenerate and non-degenerate phonons. We considered two 
forms of electron-phonon coupling: conventional, in which the phonon displacement 
couples to the electron charge density, and Jahn-Teller, in which the phonon displace- 
ment couples to the splitting of the electronic levels. We showed that as T — > three 
regimes may be distinguished. In weak coupling the mean square displacement of 
the oscillator coordinate from its non-interacting reference position vanishes and the 



27 



momentum-integrated electron spectral function assumes the non-interacting form. 
Corrections lead to the usual linear T resistivity. In intermediate coupling the mean 
square displacement of the oscillator coordinate is nonvanishing and the spectral func- 
tion is changed from its non-interacting form, but the density of states at the Fermi 
level is nonzero so the T — > resistivity is neither zero nor infinite. In strong coupling 
the mean square displacement is non-vanishing and is sufficiently large that a gap 
appears in the spectral function and the resistivity diverges as T — > 0. Examples 
of the evolution of the T — Q momentum-integrated spectral function with coupling 
strength are shown in Figs. 1 and 10. Figure 1 is calculated for the simplest model 
at half-filling, but the results at any filling for any model involving coupling of the 
oscillator to the charge density of the electrons are similar. The only important ef- 
fects are that away from half-filling in the strong coupling limit the relative weights 
of the two peaks are different and that there is at T = a discontinuous, rather than 
continuous, change in the spectral function and (r^) as the coupling is increased. We 
believe, but have not proven, that at any T > the change is continuous. The evolu- 
tion of the spectral function for the case of Jahn-Teller couphng is shown in Fig. 10. 
The significant difference is the presence of the mid-gap feature in the strong-coupling 
limit. This has relative weight |2 — 2n| and corresponds physically to sites unoccupied 
by electrons. This mid-gap feature means that for Jahn-Teller couplings and n ^ 1 
the gap controlling physical properties is much smaller at given coupling than it is in 
models in which the phonon is coupled to the charge density, and the mid-gap state 
is absent. 

Three physical effects have been left out of our calculation. One is the on-site 
Coulomb repulsion, which is certainly important in transition metal oxides. The 
Coulomb interaction tends to localize electrons, with the effect being strongest at 
commensurate fillings and so will tend to reinforce the localizing effect of the electron- 
phonon interaction. Because the localizing effect of the Coulomb interaction is 
strongest at commensurate fillings, this will also lead to an interesting doping de- 
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pendence: the effective electron-pfionon coupling will weaken as n is varied from 1. 
The Coulomb interaction will also change the form of the spectral function, because 
it raises the energy of states with two electrons on the same site. For example, the 
uppermost peak of the strong-coupling spectral function shown in Fig. 10 corresponds 
physically to the density of states for adding an electron to a site with one electron 
already present, at fixed phonon configuration. Coulomb effects will shift this energy 
upwards. A quantitative treatment of these effects within the present formalism will 
be left to a future paper. 

A second important piece of physics is quantum fluctuations of the phonons. At 
low temperatures and in the absence of long range order, these will allow electrons to 
move from site to site even in the strong coupling regime, so cutting off the strong- 
coupling divergence in the resistivity. One mathematical consequence is the appear- 
ance of Gaussian tails to the spectral functions in the gap regions. 

A third omitted piece of physics is long range order. It has been assumed through- 
out that the lattice distortions are random from site to site; thus the localization is 
due purely to polaronic effects. Near commensurate densities, charge density wave 
effects will be important as well. 
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APPENDIX A: RELATION TO MIGDAL-ELIASHBERG THEORY 

In this Appendix we derive the relationship between the formalism we have used 
and the conventional Migdal-Ehashberg theory of the weakly coupled electron-phonon 
system. The essential parameter of the conventional treatment is a dimensionless 
coupling constant Xconv which is defined in terms of the electron self energy E(a;,T') 
via 
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dE{u,T = 0) , 
hm^^o = Konv (Alj 

To establish the relationship we introduce a finite phonon mass Mph- This implies 
a Debye frequency ujd given by 

ul = k/Mpho" (A2) 

with a the lattice constant. (We use units in which h= 1.) We then introduce phonon 
creation and annihilation operators b\b via 

r = {b^ + b)/{2MphUJDY^^ (A3) 
Comparison to Eq. H yields 



(A4) 



The standard electron-phonon calculation gives 



— = norbnspin-^T^K^T) (A5) 

where 'D{ejr) is the single-spin density of states at the Fermi surface. Thus 

Aeon. = tV{er)\ (A6) 

Elsewhere in this paper it has been shown that in the dynamical mean field cal- 
culations polaron effects occur at A ~ 1 implying \conv ~ 1- For example, in the 
model with Uorb = ngpin = 1 and a semicircular density of states at /i = 0, frozen 
phonon distortions begin to occur at A = vr {\conv = !)• The critical values of A found 
here are articifically small because of neglect of quantum fiuctuations and intersite 
interactions. In models with quantum phonons, somewhat larger A values will be 
required to produce insulating behavior, but the general conclusion that metallic be- 
havior breaks down at a Xconv not too much larger than 1 will still hold. Now Migdal 
showed, using phase-space arguments, that the parameter which controls perturba- 
tion theory about the zero electron-phonon coupling limit is not A but Xuj£,/t (or if 
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Therefore, one might expect to be able to use the Migdal-Ehashberg 
equation to study the crossover from Fermi hquid to polaron physics. The difficulty 
with this argument, however, is that the ground state about which to perform the 
expansion is not known for A > 1. In fact, in the limit considered in this paper, the 
Migdal-Ehashberg equations are identical to those obtained by expanding Eq. ^ to 
order and solving self-consistently. Corrections to this self-consistent solution are 
indeed small, if XT/t is small, but the starting point is seen to be wrong if A > Ac ~ 1. 

APPENDIX B: NON-DEGENERATE FERMIONS WITH LORENTZIAN 

DENSITY OF STATES 

In this Appendix we use steepest descents techniques to analyze the T ^ limit 
of spinless fermions with a Lorentzian density of states coupled to a scalar classical 
oscillator. A similar technique was used in Re^. The partition function is 

Z{g,n) = J dre-^^''^'''^/^ (Bl) 

with 

E{r, g, n) = -r^ — grn — [/U — gr][- H — tan~^(n — gr)] + —ln[l + (/x — gr^] + jjn + OT^ . 
2 2 TT 27r 

(B2) 

and the chemical potential /i chosen so that 

dZ/dii = 0. (B3) 
As T — > 0, the r integral is dominated by the values r* minimizing E. We have 

Z{g, n)^Y. Zaexp - E{rl g, n)/T. (B4) 

a 

Here a labels the extrema; r* is a solution of 

1 1 

r = g{n tarT^ [fi — gr)) (B5) 

2 TT 
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and > is the weight associated with extremum a. Each Za is a product of two 
contributions, one from integrating over the quadratic fluctuations in r about r = r* 
and one from the leading T-dependence of In the present problem we have, at low 
T, fi = fiQ + AT; the T-linear term contributes to Za. 

Equation BSI may have one or three solutions. In the latter case, either one has 



the lowest energy or two are degenerate. If there is one dominant extremum then Eq. 
p3| implies 

fi = gr + tamrln — 1/2). (B6) 

so Eq. implies r = 0. 

Now consider the case of two degenerate extrema, at r = ri and r — r2. From Eq. 



B3| and Eq. |B5| one finds ZiVi + Z2r2 = 0; as Zi^2 > 0, ri and r2 must have opposite 
signs. It is convenient to define 

R =^(^^ + ^2)/2, 
A =g{ri-r2)/2. 
In terms of these variables Eqs. become 



(B7) 



R = /(„ - 1/2) - 'tjan'^-^t^l^^, (B8) 
. 1 2A 

A = —tan-^ — . (B9) 

271 1 + (/i - i?)2 - A2 ^ ^ 

These equations must be solved subject to the constraint E{ri,g,n) = E{r2,g,n) 



which by use of Eqs. |t38i [B2| may be written 



q"^ , 1 + (R- u + 



At small R and A, Eqs. |B1CI| and are only consistent if R = fi = g^{n — 1/2), so 



B10| is trivial. This may be most easily seen by expanding the two equations in A 
at fixed ^ — R and comparing the A^ coefficients. This immediately implies that if 
n ^ 1/2 the transition is first order. 

In the strong coupling limit, ri = —g{l — n) , r2 = gn and /i = \g'^{n — 1/2). 
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APPENDIX C: NON-DEGENERATE FERMIONS WITH SEMICIRCULAR 

DENSITY OF STATES 

In this Appendix we analyze the T — hmit of non-degenerate fermions with a 
semicircular density of states coupled to a classical oscillator. The technical differences 
from the Lorentzian case treated in Appendix B are that the mean field parameter 
satisfies the self-consistency equation Eq. |TH] and that the energy is given by 

E{r,g,n) = -|TX;n[«n - {iuJn + /u)]^ + T~ dm -TY.nln[an - gr] + fm. (CI) 
with r given by a solution of 

r = g{n + j:^^)- (C2) 

If there is one dominant extremum then Eqs. |B3| and |18| imply n = — J2n{(^n + 9^)^^ 
so Eq. implies r = 0. In this case 

a-i^^n) = hii^^n + At) - sgnun\j {iujn + -4]. (C3) 



If there are two degenerate extrema, Eqs. p^ , p7| apply, as do the relations 
ZiVi + Z2r2 = and E{ri, g, n) = E{r2, g, n) used in Appendix B. It is convenient to 
rewrite the equations in terms of i?, A, 6„ = a„ + and Zn = ioJn + /i + -R as 



^ = -fTY.j^y (C5) 



2R^ = 2g'n^ + g^TY.ln^±^, (C6) 
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For n 7^ 1/2, the transition may be shown to be first order by the argument used 
in Appendix B: one expands Eqs. and to order at fixed b and R and 
observes that the equations are not compatible. 

For n = 1/2, fi = R and b is odd in iun, so Eqs. |C^, ^ are satisfied trivially. 



The other two equations may be easily solved by viewing them as equations for b 
and at given A. Because Eq. |C7| is only cubic, the solution may be written down 
immediately. It is easist to find the proper branch of the solution if Z is continued 
to the real axis. Once b is found, may easily be computed from Eq. |C^. From 
Eq. one sees immediately that the imaginary part of 6, b" , vanishes at = for 
A^ = 1. The coupling corresponding to this, g = 1.63 is the one at which a gap first 
appears in the spectrum. Similarly, one may show that the leading (u;^) correction to 
the uj = value of b" vanishes at A^ = 1/4. 

The strong coupling limit is analytically tractable at all n. b" is nonzero only for 
Z near ±A. Defining b± = b± A and Z± = Z ±A and neglecting terms of order 1/A 
one finds from Eq. 

6± = Z±-^(l±i?/A). (C8) 

This equation may be solved and the results inserted into the other equations. One 
finds 

ri = -g{l - n), 

r2 = gn, (C9) 
/X =g\n- 1/2). 



The physical content of these solutions is one band, of width v2n centered at 
UJ = —g'^ representing occupied states and another ui = g"^ of width ^2(1 — n) 
representing unoccupied states. 
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APPENDIX D: ORBITALLY DEGENERATE ELECTRONS WITH 
LORENTZIAN DENSITY OF STATES 



In this Appendix we analyze the T ^ hmit of orbitally degenerate electrons 
with a Lorentzian density of states coupled to a classical Jahn- Teller oscillator whose 
displacement has magnitude r. The treatment parallels Appendix B but the details 
are different. Instead of B-2 we have 

E{r,g,n) = Jr^ - {fi + gr){l- + -tan~\fi + gr)) 
-(/^ - 9r){- + -tan~^{fi - gr)) + ^ln[l + {fi + gr)^] 

Z 71 Z7l 

+ ^ln[l + (/i - gr)^] + fin 

ZTT 

with n chosen so B-3 holds. Instead of Eq. ^ we have 

^- [tan~^ (/i + gr) — tan^^ (/i — gr)] . (Dl) 

Suppose first there is only one dominant extremum. If g is small this is at r = 
and 

H = tan^{n-1) (D2) 

and 

2 TT 

E{r = 0, ra) = ln{cos-{n - 1)). (D3) 

71 2 

At g = gdfi) = 2/(7r(l + /i^), there is a second order transition to a state with r ^ 0. 
In this state, the chemical potential is fixed by 

1 



r = 

TT ■ 



n 



1 + + gr) + tan~^{fi - gr)]. (D4) 



As (7 — >■ cxD at n < 1, this solution tends to r = igfn, /i = tan7r{n — 1/2) — g'^n, 
E = —\g'^n'^ + Hng'^n. 

At large g and n 7^ 1, Eq. pi| has an alternative solution, r = g. This is incompat- 
ible with Eq. p4| , so the solution must be degenerate with another extremum which, 

35 



by symmetry, must be at r = 0. These requirements imply that /i = \g^sgn{n — 1) 
and E — > —^g'^n + ^Ing. If n is near 1 the two extremum solution becomes lower in 



energy at (7 ~ 1 — n ^ g^, if n is near 0, the two extremum solution is lower in 
energy at all g > gc- 

APPENDIX E: JAHN-TELLER COUPLED ELECTRONS WITH 
SEMICIRCULAR DENSITY OF STATES 

In this Appendix we analyze the T limit of orbitally degenerate electrons 
with a semicircular density of states coupled to a classical Jahn- Teller oscillator whose 
displacement has magnitude r. The treatment follows Appendix C but with technical 
differences analogous to those found in Appendix D. Instead of Eq. |C1| we have 



E{r, g, n) = T 5][a„ - + /i)]^ + Ir^ - T ^ ln[al - gh^] + /in. (El) 



Instead of Eq. IC2| we have 



If there is one dominant extremum then a satisfies 



while /i is given by 



a„ = + /i „ „ (E3) 



^ al- g^r^ 



For small g, the only possible solution is r = 0, so a„ is given by Eq. and /i may 
be found numerically from Eq. The r = solution becomes linearly unstable at 

g'M = ^^/i^~f^'f' (E5) 

to a solution with r 7^ 0. This solution implies a frozen in Jahn- Teller splitting 
random from site to site but of the same magnitude on all sites, and therefore a two 
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peaked spectral function. Again the equations simplify in the strong coupling limit, 
where the reasoning that led to Eq. ^8| shows that the single-extremum solution 
tends to one with a spectral function consisting of two semicircular features of half 
width (vs half- width 2 at g < gc), centered at = — /i ± g'^n. For n < 1 yU is 
such that the lower band is partly filled, thus fi = —g'^n + and the energy is 

— \g'^n'^ + 0{1 — n). the (9(1 — n) term comes from carrier kinetic energy in the partly 
filled band. 

As in Appendix D, an alternate two-extremum solution exists, with one extremum 
at r = and the other at r 7^ 0. As (7 — > cxd, the nonzero r ^ g. The spectral function 
consists of three well separated peaks; two each of weight min{n, 2 — n), at u — fi = ±g 
corresponding to sites with a electron and one, of weight 2|1 — n|, centered at = — /z 
and corresponding to the empty sites if n < 1 and the doubly occupied sites at n > 1. 
The chemical potential is —^g"^ and the energy is —\g'^n + Ol/g. As (7 ^ 00 the two 
extremum solution is lower in energy. For n the divergence of gc shown in Eq 
[E]5| implies that the undistorted state proceeds directly to the three peaked one via 
a first order transition. For n near 1 the calculation sketched in the text shows that 
the two-peaked solution, with a gap between the two peaks, exists for a finite range 
of g 
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APPENDIX: FIGURE CAPTIONS 



Fig. 1 Momentum-integrated electron spectral function A{uj) plotted against frequency 
CO for non-degenerate electrons coupled to non-degenerate phonons at density 
n — 1/2, temperature T = and several different values of frozen- in lattice dis- 
tortions r — (heavy dashed hne), .87 (heavy dotted hne), 1.0 (hght solid line), 
1.17 (light dashed line) and 1.49 (light dotted hne) corresponding to = 1.53, 
1.60, 1.63, 1.69 and 1.83, respectively. Inset: coupling constant (A) dependence 
of r. 

Fig. 2 Temperature (T) dependence of resistivity p calculated for non-degenerate elec- 
trons coupled to non-degenerate phonons at density n = 1/2 for the coupling 
constants used to create Fig 1: g = 1.29 (heavy solid line), 1.53 (heavy dashed 
line), 1.60 (heavy dotted line), 1.63 (light solid line), 1.69 (light dashed line), 
and 1.83 (hght dotted line). Note that g — 1.29 and g — 1.53 have the same 
lattice distortion r = at T=0 and therefore the same spectral function. 

Fig. 3 Temperature (T) dependence of electron spectral function A{uj) for g = 1.69, 
n = 1/2 and T = (light sohd hne), .01 light dashed hne), .02 (hght dotted 
hne), .07 (heavy sohd) and .15 (heavy dashed line). Inset: spectral function 
plotted over a wider range of frequency for parameters used in main figure. The 
higher the T the longer the high frequency tail. Comparison to the appropriate 
curve in Fig. 2 shows that the rise in the resistivity begins at the T — 0.07 at 
which A first begins to develop a minimum at a; = 0. 

Fig. 4 Temperature dependence of phonon probability distribution P{r) for g — 1.69, 
n = 1/2, and T = .01 (light dashed line), .02 (hght dotted hne), .07 (heavy 
solid line) and .15 (heavy dashed line). 

Fig. 5 Low temperature behavior of resistivity for the three models considered in the 
text for parameters such that the spectral function has a gap at T = 0: lower 
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panel: non-degenerate phonons and non-degenerate electrons, middle panel: 
non-degenerate phonons and degenerate electrons, upper panel: degenerate 
phonons and degenerate electrons. The expected low-T behavior is p ~ j^x^a/t 
with A the low-T gap in the spectral function. The heavy dots mark values of 
A obtained from an analytic T = calculation as described in the text. 

Fig. 6 Temperature (T) dependence of resistivity (p) for model of degenerate electrons 
and non-degenerate phonons with n = 1 and g = 0.91 (heavy solid line), 1.09 
(heavy dashed line), 1.13 (heavy dotted hne), 1.15 (light solid hne), 1.20 (hght 
dashed line) and 1.30 (hght dotted hne). These couplings produce the same 
T = gaps as in Fig. 2. The factor of \/2 in the couplings and the factor of 2 
in the resistivity relative to Fig. 2 comes from orbital degeneracy. 

Fig. 7 Temperature (T) dependence of resistivity (p) for model of degenerate electrons 
and degenerate phonons with n — 1 and g — 0.91 (heavy solid line), 1.09 (heavy 
dashed line), 1.13 (heavy dotted hne), 1.15 (light sohd hne), g = 1.20 (light 
dashed line) and 1.30 (light dotted line). These couplings produce the same 
T — gaps as in Figs. 2, 6. Note p is much larger than that shown in fig 6 for 
the non-degenerate electrons model, reflecting the stronger scattering. 

Fig. 8 Lower panel: phonon probability distribution P(r) for g — 1.085, n = 1, T = .05 
(heavy solid line), .1 (light solid line) and .15 (light dashed line) and non- 
degenerate phonons. Upper panel: phonon probability distribution P(r) for 
same parameters and degenerate phonons. 

Fig. 9 Comparison of temperature evolution of electron spectral functions for degen- 
erate (upper panel) and non-degenerate (lower panel) phonons and orbitally 
degenerate electrons for n = 1 and g = 1.085, T = .05 (heavy solid line), .1 
(hght sohd line) and .15 (hght dashed hne) as used in Fig. 8. 

Fig. 10 Evolution of T = spectral function with coupling strength. Model with de- 
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generate electrons, degenerate phonons and n — 0.75 for coupling constants 
= 1.53 (light solid line), 1.58 (light dashed hne), and the strong coupling limit 
solution evaluated aX. g = 2.47 (heavy solid line). 

Fig. 11 Comparison of low T spectral functions. Upper panel Jahn- Teller (heavy line, 
shifted hy ^ — —1.23) and non-degenerate phonons (light line) at T = .03 
g — 1.58, and n — 0.75. Lower panel: n — 1 (light) and n — 0.75 (heavy) 
with Jahn- Teller phonons, T = 0.03 and g = 1.58. n = 0.75 curve (heavy line) 
shifted hy fjL — —1.23. 

Fig. 12 Comparison of low T phonon probability distribution P(r) for Jahn- Teller 
(heavy line) and non-degenerate (light line) phonons for T — .03, g — 1, and 
n = 0.75. Note that the large r peaks occur at almost exactly the same r. Note 
also that the Jahn- Teller P(r) is defined only for r > 0. The asymmetry in 
P(r) in the non-degenerate case is a consequence of particle hole assymmetry 
due to n 7^ 1. 

Fig. 13 Comparison of temperature dependent resistivities for Jahn- Teller (heavy hne) 
and non-degenerate phonons (light line), g = 1.29 and n = 0.75. 
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